Titanium-capped carbon chains as promising new hydrogen storage 
media 
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The capacity of Ti-capped sp carbon atomic chains for use as hydrogen storage media is studied using first-principles density 
functional theory. The Ti atom is strongly attached at one end of the carbon chains via d-p hybridization, forming stable TiC„ 
complexes. We demonstrate that the number of adsorbed H2 on Ti through Kubas interaction depends upon the chain types. For 
polyyne (n even) or cumulene (n odd) structures, each Ti atom can hold up to five or six H2 molecules, respectively. Furthermore, 
the TiCs chain effectively terminated on a C20 fullerene can store hydrogen with optimal binding of 0.52 eV/H2. Our results 
reveal a possible way to explore high-capacity hydrogen storage materials in truly one-dimensional carbon structures. 
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1 Introduction 



To develop economical hydrogen energy, carbon nanos- 
tructures with sp 2 -like bonding functionalized by transition 
metal (TM), 1-6 alkali metal (AM) 7,8 and alkaline-earth metal 
(AEM) atoms 9 have been expected to be promising storage 
materials due to their light weights and large surface areas. 
Recently, Iijima et a/.— transformed graphene to single sp- 
hybridized carbon chains containing 16 carbon atoms. In ad- 
dition, chains with metal atoms connected to the ends have 
been previously generated, 11 and their magnetic, electronic, 
and transport properties have been studied extensively. 12,13 
These results inspired us to consider whether metal-capped 
carbon atomic chains are efficient hydrogen storage media. In 
this study, we conduct theoretical studies of high-capacity hy- 
drogen storage media consisting of a Ti atom capped on two 
kinds of atomic carbon chains, cumulene (with double C-C 
bonds) or polyyne (with alternating singlet and triplet bonds). 
The number of adsorbed H2 molecules is only determined by 
the type of chain. Each Ti atom in TiCs and TiCs can bind 
up to six and five H2 molecules, respectively, corresponding 
to storage capacities of 10 wt % and 6.5 wt %. The aver- 
age binding energy of H2 molecules on TiCs (TiCs) is 0.59 
eV/H2 (0.57 eV/H2), which is between the physisorption and 
chemisorption energies. 
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2 Computational method 

Numerical calculations were performed using spin-polarized 
density functional theory (DFT) with the Perdew-Wang 
(1991) exchange-correlation function,— as implemented in 
the DMol 3 package (Accelrys Inc.).— A double numerical 
atomic orbital augmented by d-polarization functions (DNP) 
was employed as the basis set. In the self-consistent field 
calculations, the electronic-density convergence threshold was 
set to 1 x 10 -6 electron/A 3 . Geometric optimization was per- 
formed with convergence thresholds of 10 -5 Ha for the en- 
ergy, 2x 10 -3 Ha/A for the force, and 10 -4 A for the atomic 
displacements. We performed normal-mode analysis of the 
obtained structures to ensure that the structures optimized 
without any symmetry constraints were true minima of the 
potential-energy surface. 

3 Results and discussion 

We first consider the bonding energetics of a single Ti atom at 
one end of the C n (^=5-16) chains. The binding energies of Ti, 
shown in Fig. 1(a), are much larger than those of Ti on sp 2 - 
hybridized carbon materials, 2-5 and the adsorption strengths 
are thus strong enough to maintain stable Ti-C n structures at 
room temperature. Moreover, the Ti binding energies exhibit 
an interesting odd-even oscillatory behavior, consistent with 
the trend of Mulliken charge [see Fig. 1(a)]. This behav- 
ior suggests that the bonding between the Ti atom and car- 
bon chains has a significant ionic characteristic. In addition, 
the configurations of T\C n complexes exhibit significant dif- 
ferences for odd or even values of n. For instance, the C-C 
distances in TiCs, depicted in Fig .1(b), are rather uniform. In 
contrast, the bond lengths of TiCs alternate, consistent with 
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G-7T bonding [Fig. 1(c)]. Contour plots of the total electronic 
charge density confirm these bonding patterns. 




Energy (eV) Energy (eV) 

Fig. 1 (a) Mulliken charge and binding energy of TiC n as a function 
of the number of carbon atoms. The optimized bond lengths (in A) 
and the contour plot for the total electronic charge density of (b) 
I1C5 and (c) TiCg. The large and small balls represent the Ti and C 
atoms, respectively Calculated partial density of states (PDOS) of 
(d) TiCs and (e) TiCs in units of state/eV. The solid black lines and 
dotted blue lines represent the C-2p and Ti-3 d states, respectively 
Insets of (d) and (e) show the bonding orbitals (the two colors 
denote ± signs of the wave function) corresponding to states 
indicated by black arrows. 

To elucidate the nature of the bonding, the projected den- 
sity of states (DOS) shown in Fig. 1(d) illustrates the bond- 
ing orbitals of the C5 (Cs) and Ti atom result from the hy- 
bridization of the carbon chain 71 orbitals and the Ti-d orbitals. 
The C-p and Ti-d hybridization is evidently stronger in TiCs 
than in TiC5 , in accordance with the larger binding energy of 
TiCs. Since odd-numbered and even-numbered free chains 
have full and half-occupied HOMOs (highest occupied molec- 
ular orbitals), respectively, different states are responsible for 
the bonding of the Ti atom to C5 or Cs. From the isosurface 
plots of these states [Fig. 1(d)], it is clearly that the LUMO (a 
linear combination of degenerate 7l Px * and n Py * antibonding 
orbitals) of C5 hybridizes with two Ti-d orbitals to create a 
new 71-symmetry bond. However, for TiCs there is an overlap 
between one Ti-d orbital and the iz Py * antibonding state of Cs. 
As mentioned above, the bonding mechanisms for TiCs an d 
TiCs seem to be very similar to the Dewar coordination.— 

The stability of the TiCs and TiCs complexes was further 
tested by normal mode analysis, and no negative mode was 



found. However, it remains to be shown whether it is possible 
to realize linear C n -Ti chains experimentally. To address this 
issue, we calculated the minimum-energy path (MEP) for dif- 
ferent reaction processes with the nudged elastic band method 
to determine the energy barrier. 17 The image number consid- 
ered is 16 to ensure that the obtained MEP is correct. In Figs. 
2(a) and 2(b), the MEPs show that a single Ti atom is attached 
to the right end of C5 and Cs without any energy barrier. The 
other end is free and ready to be attached to any other struc- 
ture with a high surface-to- volume ratio. The main Ti vibra- 
tional frequencies, 623 cm -1 in TiCs and 359 cm -1 in TiCs, 
correspond to the stretching modes of C-Ti. These character- 
istic modes will provide a reference for the Raman/IR spectra 
of the synthesized materials. To this end, a monatomic car- 
bon chain, realized by removing the carbon atoms from the 
graphene through energetic electron irradiation inside a trans- 
mission electron microscope, 10 can be used as the initial ma- 
terial in the fabrication process. 
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Fig. 2 Variation in the total-energy difference of a Ti atom attached 
at a distance d from the right free end of the carbon chain. 

Now we investigate the interaction between these com- 
plexes and hydrogen molecules. The first H2 molecule is 
found to dissociate on both TiCs and TiCs complexes and 
bind atomically to Ti, with interatomic H distances of 3.003 
A [Fig. 3(a)] and 2.748 A [Fig. 3(b)], respectively. Such a 
dissociation is observed in most TM-hydrogen binding [Ref. 
2] due to excessive charge transfer from the TM to the anti- 
bonding a* state of H2. Additional H2 molecules, however, 
bind molecularly around the Ti atom, since the charge transfer 
per H2 molecule is not enough to destabilize the dihydrogen 
state when more hydrogen molecules are added to the system. 

For TiCs, several initial configurations of the second H2 
molecule were considered in the search for the lowest energy 
structure. The structure in which one of the H2 molecules 
binds as a dihydride is 0.09 eV higher in energy than the 
structure in which both H2 molecules bind in molecular form 
[Fig. 3(a)]. We continue to introduce successive hydrogen 
molecules near the Ti atom. Surprisingly, it is energetically 
favorable for the TiCs complex to adsorb 6 H2 molecules, cor- 
responding to a ~10 wt % gravimetric density. A close exam- 
ination of the geometry of the TiCs(H2)6 configuration shows 
that the adsorbed hydrogen molecules are divided into two sets 
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[Fig. 2(a)]. First, the molecule on the top binding site, which 
is vertical to the carbon chain axis, has a bond length of 0.808 
A and the largest bonding energy, 0.82 eV. Second, the H-H 
bonds of the side H2 molecules lie almost parallel to the chain 
axis, with the lower H atoms of each H2 molecule tending to 
tilt toward the chain. The bond lengths between the Ti and 
the lower (upper) H atoms are 1.910 A (1.970 A). The side 
molecules are more weakly bound, with a binding energy of 
0.54 eV/H 2 and an elongated H-H bond distance of 0.810 A, 
on average. All of these features suggest that different bond- 
ing mechanisms (to be discussed below) for the two sets of H2 
molecules on TiC5 . 

On the other hand, when the second H2 molecule is bound 
to TiCg, the structural configuration shown in Fig. 3(b) is more 
energetically favorable than its isomer (not shown) by 1 .02 eV. 
As we add a third H2 molecule close to the Ti atom, all of the 
hydrogen molecules are almost parallel to the chain and are 
bonded molecularly. Interestingly, the H-H bonds turn nearly 
perpendicular to the chain when the fourth H2 molecule is at- 
tached. The maximum number of H2 molecules that could be 
bound to the TiCs complex is 5, which corresponds to a hy- 
drogen wt % of 6.5. Note that the adsorption of the fifth H2 
molecule does not affect the degree of the bond-length elon- 
gation of the side hydrogen molecules [Fig. 3(b)]. The H-H 
bond of the top H2 molecule is 0.756 A and the bond length 
between this H2 molecule and the Ti atom is 2.334 A, leading 
to an adsorption energy of 0.29 eV, which is only about half 
the value of the side molecules (0.66 eV). 

Insight into the nature of Ti-H2 bonding and the orienta- 
tion of H2 molecules can be gained from the DOS and molec- 
ular orbital (MO) analyses. The metal-dihydrogen binding 
of TiC5(H2)6, shown in Fig. 3(c), stems primarily from the 
Kubas interaction. 1 - In the energy range from -10 to -6 eV 
(region I), the d orbitals of Ti are hybridized with the a or- 
bitals of the hydrogen molecules, resulting in charge transfer 
from the hydrogen to the metal. Region II of Fig. 3(c) high- 
lights that the hybridization of the Ti-d orbitals with the H2 
G*-antibonding orbitals is responsible for keeping the side H2 
molecules parallel and the top H2 molecule perpendicular to 
the TiC5 chain. The isosurface plots of these states, as shown 
in Fig. 3(d), clearly confirm that the bonding is dominated 
by the overlap between the Ti-d and H2-S states. To better 
understand the orientation of the side H2 molecules, we also 
plotted the deformation charge in Fig. 3(a). The charge distri- 
bution around the lower H sites, marked in yellow, indicates 
that these H atoms lose a greater amount of charge. With the 
help of Mulliken charge population analysis, we identified that 
the upper H atom has a charge of 0.04 e, while the lower one 
carries 0.12 e due to charge transfer from H to the C atom. 
Thus, the distances between the Ti and the upper or lower H 
depend on the strength of the Coulomb interaction between 
the H and C atoms. The Kubas interaction is likewise found 
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Fig. 3 Optimized configurations of H2 molecules adsorbed to the Ti 
atom on (a) C5 and (b) C% along with the typical bond lengths (in 
A). The large, medium and small balls represent the Ti, C and H 
atoms, respectively. The right panels of (a) and (b) show the 
deformation electron densities (molecular charge densities minus 
atomic charge densities) for the TiC5(H2)6 and TiCs (112)5 
complexes. The deformed density marked in blue corresponds to the 
region that contains excess electrons, while that marked in yellow 
indicates electron loss. The isovalue equals 0.1 e/A 3 . (c) PDOS for 
six hydrogen molecules attached to the Ti atom on C5 . (d) 
Isosurfaces of major molecular orbitals corresponding to states in 
region I and II of panel (c). Panels (e) and (f) show the a and a* 
orbitals, respectively, for the side hydrogen molecules hybridized 
with the Ti-d orbitals in TiCs 0^2)5. (g) The bonding orbital for the 
top hydrogen molecule. 



in TiCg(H2)5. Figure 3(e) shows that the Ti-d orbitals inter- 
act with the a orbitals of the H2 molecules. The hybridiza- 
tion between the Ti-d orbitals and the side H2 molecules' a*- 
antibonding orbitals plays a role in keeping them planar [Fig. 
3(f)], in contrast to the situation in TiCs(H2)6. The top hy- 
drogen in TiCg(H2)5 interacts with Ti rather weakly compared 
to the rest of the hydrogen molecules, as the relevant MO is 
more heavily polarized toward the Cs than toward the H2 [see 
Fig. 3(g)]. The Kubas interaction mentioned above leads to 
variations in the effective charge of the Ti atom. The Mul- 
liken charge of the Ti atom varies from positive to negative as 
the number of hydrogen molecules increases. When the fifth 
(sixth) H2 is adsorbed onto TiCs (TiCs), the Ti atom carries 
-0.13 (-0.49) electrons. 

The above results demonstrate that the interaction between 
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one Ti atom and the carbon chain is very strong, and the bind- 
ing energy of H2 on TiC n is optimal for room temperature ap- 
plication. Now the question arises: can the storage perfor- 
mance of TiC,2 be influenced by introducing another Ti? We 
first discuss the second Ti atom saturating the free end of TiC n 
(n= 5, 8). The binding energy for Ti to the TiCs (TiCs) is 4.31 
eV (4.39 eV), which is obviously much large than that of Ti on 
the sp 2 carbon nanostructures. 2-5 However, further function- 
alization through Ti atoms bridging other C-C bonds makes 
the chain unstable, due to the fact that the small separation be- 
tween Ti atoms leads to stronger Ti-Ti coupling. At the end, 
the linear chains are transformed to the clusters composed of 
C and Ti atoms, which is also confirmed by other theoretical 
studies. 19 We therefore only consider two Ti atoms each cap- 
ping one end of the C n . 

Next, we studied the H2 storage capacity of the TiC^Ti 
complex. The TiCsTi and TiCsTi can adsorb 12 and 10 H2 
molecules, respectively, with average binding energies of 0.53 
and 0.55 eV/H2. As shown in Figs. 4(a) and 4(b), the gravi- 
metric densities of stored H2 molecules for TiC^Ti becomes 
13.3 and 9.4 wt % when n=5 and n=8, respectively. Analysis 
of the electronic structures of the TiC5Ti-(H2)i2 and TiCsTi- 
(Fbho shows that the Kubas interaction between H2 and TiC„ 
mentioned above also applies to the TiC^Ti complex. Note 
that introduction of the second Ti atom to TiC n not only keeps 
the number of adsorbed H2 per Ti on odd- or even-numbered 
chains, but also remains the binding energies of the H2 with 
metal atoms . 
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of adsorbed H2 molecules is six and five, respectively. The 
corresponding bond lengths of H2 and Ti-H2, as depicted in 
Fig. 5(a) and 5(b), are very similar to those in TiCs(H2)6 and 
TiCs(H2)5. Furthermore, the average binding energies of the 
H2 molecules to TiQs and TiCi6 are, respectively, 0.54 and 
0.59 eV/H2, which are very close to the values found in the 
case of TiC5(H2)6 and TiCs (£[2)5. The bonding mechanisms 
between Ti and H2 presented above also hold for these chains. 
This suggests that the length of the chains does not affect the 
hydrogen adsorption performance. In summary, the finding 
that a single Ti atom-capped, odd-numbered ( even-numbered) 
carbon chain can bind up to six (five) hydrogen molecules is 
very general. 
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Fig. 5 Geometries of (a) TiC 15 (H 2 ) 6 , (b) TiC 16 (H 2 ) 5 , and (c) a 
TiCs chain terminated on the C20 fullerene. (d) Optimized geometry 
of 6 H2 molecules adsorbed on C2o(TiC5). The large, medium and 
small balls represent the Ti, C and H atoms, respectively. 
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Fig. 4 Geometries of (a) TiC 5 Ti-(H 2 )i 2 and (b) TiC 8 Ti-(H 2 )i . The 
large, medium and small balls represent the Ti, C and H atoms, 
respectively. 

Because controlling the number of atoms per chain in the 
synthesis of carbon chains may be difficult with today's tech- 
nology, it is important to know whether the results reported 
above for TiCs and TiCs are applicable to other Ti-capped car- 
bon chains and how they vary with chain length. To this end, 
we have also studied the potential of the TiQs and the TiCi6 
complexes as storage media, and we found that the number 



While the above results are promising for isolated TiC„ 
systems, one can imagine terminating the other end of these 
complexes with suitable graphitic nano-objects. These new 
structures not only represent the typical interface in realistic 
nanostructures produced by cluster beam deposition. 20 Here 
we choose C21 fullerene as the end-capping candidate, as it 
is the most frequently experimentally synthesized curved sp 2 
system. An increasing amount of experimental 21 evidence 
shows that these hybrid sp+sp 2 carbon-based systems (with 
linear chains connecting s/? 2 -type fragments) exhibit unusual 
electronic and optical properties. Figure 5(c) illustrates that 
a TiCs chain can be effectively stabilized by termination on 
C20. The calculated binding energy (with respect to the iso- 
lated TiCs plus the fully relaxed C20) is 2.79 eV. As shown in 
Fig. 5(d), the TiCs assembled on the C20 structure can hold 6 
H2 molecules, with an average binding energy of 0.52 eV/H2. 
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4 Conclusion 

In conclusion, using all-electron DFT calculations, we have 
shown that each Ti atom adsorbed on even or odd-numbered 
carbon atomic chains can bind up to five or six hydrogen 
molecules, respectively. Note that the number of adsorbed H2 
molecules depends only on the type of chain. We propose that 
the T1C5 chain terminated effectively on a C20 fullerene can 
also store 6 H2 molecules with an average binding energy of 
0.52 eV/H2. Recent experiments have produced junctions be- 
tween a single carbon chain and two fullerenes, 22 which has 
provided a way to synthesis the sp+sp 2 system we proposed 
here. We hope that the theoretical results presented here will 
provide a useful reference for the design of high-capacity hy- 
drogen storage materials in the laboratory. 
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